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Abstract 

We speculate that stationary force-free pulsar magnetospheres are screw unstable, and the spin- 
down power of real pulsars is carried by turbulent Poynting jets. The turbulent jet entrains poloidal 
flux as it propagates away from the star. Due to poloidal flux entrainment, poloidal field inside the 
jet decreases at a slower rate than the dipole's r -3 , and the pulsar power increases accordingly: 
L ~ (/i 2 $7 4 /c 3 )(f2-R/c)~ a , where u is the magnetic dipole, is the frequency, R is the neutron star 
radius, and a < 1 is an index which depends on spin-dipole angle. Our speculation is of interest 
because it seems to provide the only possible explanation (without fine-tuning) of observed pulsar 
braking indices. 
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I. INTRODUCTION 



Pulsars spin down. The spin-down power is probably carried by Poynting jets propagating 
along open magnetic field lines, and the resulting spin-down power is thought to be 

L = fi«f-f- (1) 

where f{9) ~ 1 is an unknown dimensionless coefficient depending on the spin-dipole angle 
9 Q]. The power is comparable to the magnetic dipole power j^J because one assumes that 
pulsars simply eject the characteristic dipole field Bi c ~ /V r ?c a t the light cylinder radius 
ric = c/Q, giving the power L ~ rf c cBf c ~ fi 2 Q 4 /c 3 . 

The conventional pulsar power formula (Q) has been confirmed recently, when the axisym- 
metric pulsar magnetosphere was calculated, giving the dimensionless coefficient for aligned 
rotators 0,0]: /(0) = 1 ±0.1. Assuming that non-zero 9 solutions are similar to the aligned 
case, it would seem that the large-scale electromagnetic background, on which all the pulsar 
emission is occurring, has basically been established. 

One, however, can argue differently. Stationary pulsar magnetospheres can be unstable. 
Then real pulsar magnetospheres are turbulent, and to account for the effects of turbulence, 
one has to modify the pulsar power formula. Within the framework of FFE (force-free 



electrodynamics, see eg. 



Q 



and references therein), the only dimensionless parameter, which 
one can use to modify the pulsar power, is QR/c << 1, where R ~ 10km is the neutron star 
radius. Since the background is self-similar, the only reasonable modification of the pulsar 
power is the power law. Thus, we will postulate that the correct pulsar power is given by 

where the dimensionless coefficient / ~ 1 and index a < 1 depend on the spin-dipole angle 
9. The negative power in equation (0) corresponds to poloidal field entrainment into the 
Poynting jet. This leads to an increased value of the field within the jet, and thus to a larger 
power. 

If one uses the anomalous pulsar power formula (j2J) to infer the surface magnetic fields of 
neutron stars, the predicted values of the fields will be lower than the standard ones. Also, 
the turbulent magnetosphere model might be needed to understand the pulsar emission. 
The most important reason to take the turbulent model seriously is because it provides the 
only viable explanation of pulsar timing observations, §2. 



II. PULSAR TIMING AND ANOMALOUS PULSAR POWER 



Pulsar spin-down may be quantified by braking indices 



n 2 n 



(3) 




m = 



ft 3 



The values of the first braking index n are less than 3 in all five pulsars where it has been 
measured Jsj]. Even more important for our purposes is the fact that for Crab (n = 2.519, 
m = 10.23 ± 0.03 from Table 1 of j^) the second braking index is m = n(2n — 1) to 
2cr. Unfortunately, this result can possibly be "contaminated by timing noise and frequent 
glitches" Q. 

What are the theoretical expectations for the braking indices Q| ? 

Anomalous pulsar power (j2J), with the assumption of no change in the spin-dipole angle 
and no change in the dipole moment, gives the spin-down equation 



in agreement with both 3 > n = 3 — a and m = n(2n — 1). 

Conventional pulsar power Q gives n = 3 for constant spin-dipole angle and dipole 
moment. One needs to assume that the angle or/and the dipole moment change. Let us 
consider these possibilities. 

The model postulating variable dipole moment requires fine-tuning, because there is no 
reason to expect that the time scale for the change of the magnetic field and the time scale 
for the spin-down should be the same. If one accepts Crab's m = n(2n — 1), the fine-tuning 
problem becomes much more severe. 

The model postulating evolution of the spin-dipole angle might seem to require no fine- 
tuning. It is known 8] that magnetic dipole rotating in vacuum aligns in about the spin- 
down time. However, real stars are not spherical tops, and real magnetic fields are not ideal 
dipoles. Due to magnetic multipoles in the torque formula (and also magnetic and elastic 
deformations of the ellipsoid of inertia of the star [9|), the time scale for spin-dipole angle 
evolution is expected to be shorter than the spin-down time by a factor ~ QR/c <C 1 (or 
even smaller factor coming from the ellipsoid of inertia). Thus, this theory cannot work 
without fine-tuning. With Crab's m = n(2n — 1), the fine-tuning problem becomes much 
more severe. 
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We are aware of only one model which requires no fine-tuning in explaining n^3. The 
model postulates that during spin-down, the outward-moving superfluid vortex tubes entrain 
the superconducting magnetic flux tubes n|. This automatically gives the right time scale 
for the magnetic field evolution. However: (i) this theory still requires fine-tuning if one 
accepts Crab's m = n(2n — 1), (ii) this theory allows both n < 3 and n > 3, depending on 
the initial field configuration. 

Anomalous pulsar power appears to be the most natural explanation of pulsar timing 
observations. To confirm or refute this model, one needs to perform a time- dependent 
three-dimensional FFE simulation of a pulsar magnetosphere. The analytical methods only 



say that screw is the most promising instability [LjJ . But even the problem of FFE stability 
of an axisymmetric pulsar, for which the precise shape of the stationary magnetosphere is 
known, looks too difficult for analytical methods (the problem is not self-adjoint; there are 
singular current layers; the equatorial current layer is not an FFE current layer). 
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